Comparison between numerical-relativity and post-Newtonian waveforms from 
spinning binaries: the orbital hang-up case 
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We compare results from numerical simulations of spinning binaries in the "orbital hangup" case, 
where the binary completes at least nine orbits before merger, with post-Newtonian results using the 
approximants TaylorTl, T4 and Et. We find that, over the ten cycles before the gravitational-wave 
frequency reaches Muj — 0.1, the accumulated phase disagreement between NR and 2.5PN results 
is less than three radians, and is less than 2.5 radians when using 3.5PN results. The amplitude 
disagreement between NR and restricted PN results increases with the black holes' spin, from about 
6% in the equal-mass case to 12% when the black holes' spins are Si/M^ = 0.85. Finally our 
results suggest that the merger waveform will play an important role in estimating the spin from 
such inspiral waveforms. 

PACS numbers: 04.20.Ex, 04.25. Dm, 04.30.Db, 95.30.Sf 



I. INTRODUCTION 

Direct observation of gravitational waves may be at 
hand. Several ground-based gravitational-wave (GW) 
detectors have reached design sensitivity [H, H, Q , and 
the recent completion of the LIGO S5 science run, where 
science data have been taken for almost two years, marks 
an important milestone in the field. With the joining of 
the GEO600 and Virgo detectors for part of this run, 
an important step has also been taken in establishing a 
global network of detectors, and data from all these de- 
tectors are currently being analyzed. 

For black-hole binary systems, which are the strongest 
expected GW sources, current GW template banks use 
waveforms calculated by analytic post-Newtonian (PN) 
approximation methods to model waves from a binary's 
slow inspiral, and perturbation-theory techniques to 
model the ringdown of the final merged black hole. It 
is not clear how well these models describe the cru- 
cial merger phase, which will produce the strongest sig- 
nal. Recent breakthroughs in numerical relativity (NR) 
01 B@ have made it possible to describe the late inspiral, 
merger and ringdown in full general relativity, and future 
GW searches will be able to use "hybrid" waveforms pro- 
duced by combining PN and NR waveforms 0, H, H, ■ 

The construction of hybrid waveforms requires an un- 
derstanding of the frequency range in which PN and NR 
waveforms overlap, and how close to merger the PN wave- 
forms cease to be sufficiently accurate. Studies so far 
have focused on nonspinning binaries 0, IHl [13, IH, 13 ■ 
The broad conclusion from that work was that up to 
a few orbits before merger, standard PN approximants 
predict the phase evolution of the binary to within an 
accuracy of about 1 radian [13, [H, [13] , the lowest-order 
PN contribution to the wave amplitude disagrees with 
numerical results by about 6% over the frequency range 
considered in numerical simulations, and the error from 
higher-order PN amplitude corrections is reduced to only 



a few percent up to 10-15 orbits before merger [isl. [l3|. 

We now consider binaries in which the individual black 
holes are spinning. Spinning binaries have already been 
simulated in manv contextsjflljlljll [H, [H, [M [Ml, 
[13, m, [H m, [M [13, m, [lalMllirbut in this work 
we present the first long (> 15 GW cycles) simulations of 
spinning binaries, and compare with PN results during 
the inspiral. As a first example from the large parameter 
space, we consider equal-mass binaries with equal spins 
that are parallel to the orbital angular momentum of the 
binary. In addition, wc compare the spinning waveforms 
with their nonspinning counterparts, and begin to ad- 
dress the question of how well the spin can be estimated 
from an observation of one of these waveforms. We find 
that, during the slow inspiral stage, spinning and non- 
spinning binaries are difhcult to distinguish; while during 
the merger and ringdown the waveforms can be clearly 
distinguished. 

We summarize our numerical methods in Section [TTl 
and also the main features of the numerical simulations 
that we performed, and in Section IIIII summarize three 
PN calculations of a spinning binary's phase evolution, 
the Taylor Tl, T4 and Et methods. In Section HVl we 
compare the NR phase evolution and amplitude from the 
(Z = 2, m = 2) mode of with the PN phase calculated 
from the TaylorTl, T4 and Et approximants, and with 
the PN amplitude at restricted (quadrupole) order. We 
also compare the numerical waveforms with each other in 
Section |V] to get an estimate of how much nonspinning 
and spinning waveforms differ in the inspiral and merger 
phases. 



II. NUMERICAL METHODS AND RESULTS 

We performed numerical simulations with the BAM 
code ^2, 33, 34]. The code starts with black- hole bi- 
nary puncture initial data [H, [s^ generated using a 
pseudo-spectral code [s^, and evolves them with the x- 
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TABLE I: Summary of grid setup for numerical simulations. 
The grid parameters follow the notation introduced in [33 |: 
see text, hmin denotes the resolution on the finest level and 
hmax the resolution on the coarsest level. The outer boundary 
of the computational domain is at approximately rmax- The 
simulations with spin S/Mf = 0.75, 0.85 use one extra level 
of mesh refinement. 



Run 


hmin 


hmax 




S/M^^ = 0, 0.25, 0.5 simulations 


X^=2[5 X 64 : 5 X 128 : 6] 
X„=2[5 X 72 : 5 X 144 : 6] 
X„=2[5 X 80 : 5 X 160 : 6] 


A//42.7 
A//48.0 
A///53.3 


12M 
32/3M 
48/5M 


774Af 
773M 
773M 


S/Mi = 0.75, 0.85 simulations 


Xr,=2[6 X 64 : 5 X 128 : 6] 
X„=2[6 X 72 : 5 X 144 : 6] 
X^=2[6 X 80 : 5 X 160 : 6] 


A^/85.3 
A^/96.0 
Af/106.7 


12M 
32 /3M 
48/5M 


774Af 
773M 
773M 



variant of the moving-puncture [38|, |39| version of the 
BSSN [13, Eii formulation of the 3+1 Einstein evolu- 
tion equations. Spatial finite-difference derivatives are 
sixth-order accurate in the bulk 3j| , Kreiss-Oliger dissi- 
pation terms converge at fifth order, and a fourth-order 
Runge-Kutta algorithm is used for time evolution. The 
gravitational waves emitted by the binary are calculated 
from the Newman-Penrose scalar 5*4, and the details of 
our implementation of this procedure are given in [32| . 

The new simulations we performed for this analysis 
are summarized in Tables U and [TTl AH simulations are of 
equal-mass binaries; the black-hole punctures are placed 
on the y-axis ai y = ±_D/2, and given momenta in the x 
direction of Px = TP- The black holes each have the same 
spin angular momentum Si oriented parallel to the total 
orbital angular momentum, i.e., in the positive z direc- 
tion. The grid setup is described following the notation 
introduced in For example, x,,=2[5 x 64 : 5 x 128 : 6] 
indicates that the simulation used the x variant of the 
moving-puncture method, five nested mesh-refinement 
boxes with a base value of 64"^ points surround each black 
hole, and five nested boxes with 128'^ points surround the 
entire system, and there are six mesh-refinement buffer 
points. As summarized in Table |TT1 simulations were per- 
formed with spins S/M'i = 0.25, 0.5, 0.75, 0.85. We also 
make use of the results from nonspinning binaries, i.e., 
Si/M"^ = 0, as reported in [ll|. 

The physical parameters are given in Table |TT1 As 
the spin is increased, the mass parameter rrii decreases. 
This is partly because this quantity parametrizes the 
mass associated with the area of the apparent horizon, 
Mah — \/ A/IQ-K, where A is the area of the horizon, 
while the total mass of the black hole is estimated by a 
variant of the Christodoulou formula l42l. 
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AHA) 



(1) 



If we are to keep Mi constant as Si is increased, then 
Mah must decrease, and therefore decreases. The de- 



TABLE II: Physical parameters for the moving-puncture sim- 
ulations: the coordinate separation, D/M, the mass param- 
eters in the puncture data construction, rrn/M, and the mo- 
menta px/M. The punctures are placed on the y-axis, and 
for all simulations the total initial black-hole mass is Af = 1. 



Simulation 




D/M 


rrii/M 


pJM 


S25 


0.25 


12.0 


0.47579 


0.083813 


S50 


0.50 


11.0 


0.43277 


0.087415 


S75 


0.75 


10.0 


0.33608 


0.091435 


S85 


0.85 


10.0 


0.25628 


0.090857 



creasing value of rrii is also due to the extra "junk" energy 
that the spin adds to the initial data: the Bowen-York 
extrinsic curvature used in the initial-data construction 
contains unwanted gravitational radiation that increases 
in amplitude as the spin is increased. That radiation adds 
to the mass of the black hole, and can be compensated 
by further lowering rrii, but only up to a point: even- 
tually there is too much junk radiation, and we reach a 
limit in the Si/Mf that we can obtain. This limit has 
been found experimentally to be about Si/Mf « 0.928 
fill HI]. A form of puncture data that permits higher 
spins has been suggested [45|, but for the present study 
we consider spins no larger than Si/Mf = 0.85, which is 
well below the limit for Bowen-York data. 

The initial momenta for quasi-circular inspiral were 
calculated using a 2.5-PN-accurate procedure based on 
the results in [4y| and outlined in [2^. We expect these 
parameters to lead to inspiral with a small eccentricity, 
and this is indeed what we see in the numerical data; the 
eccentricity is typically of the order of e « 0.006. We 
found that the procedure described in [i^l to produce 
lower-eccentricity inspiral, although applicable for spin- 
ning binaries, does not yet include sufficient accuracy in 
the spin terms to be of use for the scenarios described in 
this work. 

We Richardson extrapolate our data with respect to 
numerical resolution and radiation extraction radius as 
described in [l^. We first split the = 2, m = 2) mode 
of the waveform r4'4 into amplitude and phase according 
to 



r*4 = A(4>{t))e-''l'^'\ 



(2) 



The amplitude is in turn written as a function of phase, 
yl(0), and this function is Richardson extrapolated with 
respect to numerical resolution. Five such functions 
are produced, one for each extraction radius. Rex — 
{50, 60, 70, 80, 90}A/. We then note that the dependence 
of the amplitude on extraction radius is modeled well by 



O 



1 

W 



(3) 



and applying a curve fit of the form ([3]) to we 
estimate A{(l)Rex —> oo). Although a method has been 
suggested to also extrapolate the phase to infinite extrac- 
tion radius in this work we again follow the proce- 
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dure described in and use the phase from the largest 
extraction radius, R^x = 90M. 

In the equal-mass case studied in [l^, we found that 
the numerical results were cleanly sixth-order convergent, 
and were therefore able to both remove the sixth-order 
error term by Richardson extrapolation, and to estimate 
the uncertainty due to higher-order error terms. The 
present simulations of spinning binaries do not show such 
clean convergence. All of the simulations exhibit conver- 
gence between fifth and sixth order. Since the code con- 
tains fourth, fifth and sixth-order elements, it is not ob- 
vious at which resolutions each error term will dominate, 
and we make the most conservative choice of using the 
highest-resolution results and estimating an uncertainty 
by assuming only fourth-order convergence. Even with 
this conservative estimate of the discretization error, we 
find that, as in the equal-mass case, the errors are any- 
way dominated by the finite extraction radii. In general 
we estimate the uncertainty in our waveform amplitude 
as less than 3%, and in the phase our uncertainty is 0.25 
radians over the frequency range that we will consider for 
comparison with PN results. 

These uncertainty estimates do not take into account 
the effect of eccentricity. In the nonspinning case stud- 
ied in [13, El], the eccentricity was e < 0.0016, and the 
error in the phase evolution due to the eccentricity was 
estimated as being well below the finite-difference and 
finite-extractdion-radii errors. Figure 15 of [l3| shows 
that the accumulated phase error of a simulation with 
e ~ 0.008 is around 0.2 radians, which is comparable to 
the numerical phase error. In the simulations presented 
in this work, the eccentricity can be as high at e ~ 0.006, 
and so significant eccentricity-induced phase errors may 
be expected. Based on Figure 15 in [l^, we estimate 
that such errors are no larger than 0.2 radians. This is a 
systematic error: the binary merges sooner if the eccen- 
tricity is increased, and therefore the accumulated phase 
disagreement A(/)e = 4ie=o — 0e>o will always be posi- 
tive. This systematic uncertainty should be taken into 
account in the comparisons we perform with PN approx- 
imants in Section IIVI The effect of the eccentricity on 
the wave amplitude is to produce oscillations in the am- 
plitude. For the simulations with higher spin, the eccen- 
tricity is larger, and the resulting oscillations are larger. 
This is clear in Figure [51 where the oscillations are visi- 
ble in the highest-spin case shown {Si/Mf = 0.75), and 
the highest-spin results are not shown because the oscil- 
lations reduce the clarity of the figure. We emphasize 
that the errors due to eccentricity may appear from our 
results to be larger for systems with larger spin, but this 
is not necessarily the case: the errors are simply larger 
when the eccentricity is large, and in the particular sim- 
ulations we have done the larger-spin configurations also 
have larger eccentricity, and therefore larger eccentricity- 
induced errors. 

Some of the features of the simulations are summarized 
in Table IIIII the mass and spin of the final black hole, 
and two quantities that allow comparison between simu- 



TABLE III: Selected global features of the simulations: the 
mass and spin of the final black hole, Mflnai and Sflnai/A/lnab 
the time AIa (in units of A'l) for the GW to evolve from 
Muj = 0.06 to its maximum amplitude; and the number of 
GW cycles ANgw between Mlo — 0.06 and the amplitude 
maximum. 



Simulation 


Mflnai 


Sflnal/A/flnj^i 


AtA{M) 


ANgw 


D12 


0.950 


0.680 


719 


11.0 


S25 


0.943 


0.757 


819 


12.7 


S50 


0.932 


0.826 


917 


14.7 


S75 


0.920 


0.896 


1040 


16.8 


S85 


0.911 


0.918 


1096 


18.0 



lations, and which demonstrate the orbital hangup effect 
first observed in numerical simulations in [l^. Starting 
from a GW frequency of Moj = 0.06, we list the time it 
takes the GWs to reach their maximum amplitude, Atj^, 
and the number of GW cycles up to the maximum. (Di- 
viding this number by two gives the corresponding num- 
ber of the binary's orbits.) The numbers clearly show 
that the binary's merger is "hung up" by the presence of 
spin. 

The mass of the final black hole was estimated by sub- 
tracting the radiated energy, as measured at the largest 
extraction radius Rex = 90M, from the ADM mass of 
the spacetime, as calculated on the initial timeslice. The 
spin of the final black hole was estimated by comparing 
a fit of the ringdown with analytic quasinormal mode re- 
sults. The values for the mass and spin are consistent 
with those already given in the literature [H, [l^, [lO] . 

Figure [1] shows the accumulated phase error between 
the nonspinning binary simulations presented in , and 
the various spin cases, for the ten cycles before the given 
spin waveform reaches a GW frequency of Mu! = 0.1. 
The GW phases are aligned such that the phase differ- 
ence is zero when Mcu = 0.1 and we also relabel the 
time so that t — a,t this point; this procedure is de- 
scribed in more detail in Section IIV Al We see that for 
Si/Mf = 0.25, the accumulated phase difference is about 
1.6 radians, while for Si/Mf — 0.85, the accumulated 
phase difference is almost five radians. These differences 
should be borne in mind when we compare with the post- 
Newtonian phase predictions. 



III. VARIOUS PRESCRIPTIONS FOR GW 
PHASE EVOLUTIONS IN PN RELATIVITY 

In this section., we provide formulae, extracted from 
Refs. gi, H IM [5l|, HI m, H, m , required to construct 
TaylorTl, T4 and Et templates for inspiraling equal-mass 
binary black holes, having their spins Si and S2 aligned 
with the Newtonian angular momentum unit vector I. 
(We retain G and c in the formulas in this section, al- 
though for the remainder of the paper we will adopt ge- 
ometrized units G = c = 1.) The aim is to provide GW 
phase evolutions under TaylorTl, T4 and Et prescrip- 
tions that include all the non-spinning contributions to 
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construction of fully 2.5PN accurate TaylorTl, T4 and 
Et templates are really applicable for spinning black hole 
binaries. 

We define spins as Si — G rnf Xi ^ii'^ = 1>2 and in 
this paper, we impose the restrictions, Si ■ I ~ +1 and 
Si -82 = +1. For compact binaries having aligned spin 
configurations, it is also possible to employ the restricted 
PN waveforms 



h{t) cx X [tf^ cos2 0(i). 



(4) 



FIG. 1: The accumulated phase difference between nonspin- 
ning and spinning binaries for the ten cycles before Muj = 0.1, 
aligned such that the phase disagreement is zero when Mui = 
0.1. The black-hole spins are, in order of increasing magni- 
tude of the phase disagreement, Si/Mf = 0.25, 0.5, 0.75, 0.85. 
Note that the differences are in some cases comparable to the 
phase disagreements with PN results that we will present in 
Section HVl 



To obtain the GW phase evolution (f>{t) in the TaylorTl 
approximant, one numerically solves the following two 
differential equations: 



d(l>{t) 
dt 



dx{t) 
dt 



- ^ 0.3/2 

GM ' 

-Mi 



(5a) 
(5b) 



3.5PN order and the spin effects to 2.5PN order. We re- 
quire, as in the case of non-spinning compact binaries, the 
dynamical (orbital) energy £ (a:), expressed as a PN series 

in terms of x = [G AdcJb/c^Y^^ , ujb{t) being the binary's 
orbital angular frequency, and the GW energy luminosity 
C{x). The fact that we are dealing with compact bina- 
ries having aligned spins implies that to the highest PN 
order considered in this paper, the spin vectors Si and 
^2 have their directions and magnitudes fixed. 

The 3.5PN accurate C{x) and 3PN accurate £{x) asso- 
ciated with non-spinning comparable mass compact bi- 
naries are extractable from Refs. [4l,[53,|5l|. The lowest- 
order General Relativistic spin-orbit coupling, appearing 
at 1.5PN order for maximal Kerr black holes, and spin- 
spin effects entering at 2PN order make contributions to 
C{x) and ^(a;) at relative 1.5PN and 2PN orders 
In the usual terminology adapted in the PN literature, 
this implies that spin-orbit coupling provides corrections 
to E{x) in terms of a;'^/^ with respect to its Newtonian 
counterpart and a similar rule applies for C{x). The 
monopole-quadrupole interactions affect C{x) and S{x) 
at relative 2PN order [52]. The next to leading order 
General Relativistic spin-orbit coupling aS^ects C{x) and 
£{x) at relative 2.5PN order [H^, E^. However, for in- 
spiraling spinning black holes, there exist corrections to 
C{x), also appearing at the relative 2.5PN order, due 
to the energy ffowing in to the black holes. For com- 
parable mass spinning black holes, the above-mentioned 
contributions are derived in Ref. 53 
tations, available in Refs. 149, W3, 5C 



For PN com 

mi, HI, m, i 



non-spinning black hole is modeled as a point-particle 
(spherically symmetric mass distribution), and a spin- 
ning black hole is approximately treated as a spinning 
spherically symmetric mass distribution. It is important 
to note that properties of black holes, like presence of 
black hole horizons, play no role in these PN computa- 
tions. Due to inclusion of black hole absorption effects, 
available in Ref. [5^, it is reasonable to argue that our 



Therefore, to obtain the GW phase evolution, relevant 
for equal mass black hole binaries with aligned spins, 
under the TaylorTl approximant that includes all the 
3.5PN accurate non-spinning and 2.5PN accurate spin 
effects, we require the following expressions for C{x) and 
£{x). 



C{x) 
'C(x)|ns 



C(x) 



£{x) 
f (a;)|NS 



'C(x)|ns + C{x)\s , 



2 0" 

Jg 

767 

71 

42 

1712 
T05" 



{ 1 



373 



- Ait X 



3/2 



.5/2, 



84 

18608019757 



59 
567' 



209563200 



355 
"64" 



In (4 V^) 



16655 
6048 



r.7/2 



(6a) 



1712 
T05"^ 

(6b) 



= 5G^^-(^^ 
+ T^XiX2 + 



1 



X2)a;3/2-H- (xi^+X2^ 



16 
1 

32 



799 
504 



(Xi +X2) 



Xi+X2+3xi +3x2^ k 



^5/2 



'S^(2;)|nS +^(2;)Is , 



Mc2 



8 

1427365 
331776 
A/c2 



37 



1069 



205 
384 



(6c) 
(6d) 



(6e) 



2 2 

X 



335 
144 



f 7 1 

|g (Xi+X2)a;'''/'-^ (X1+X2) 

{xi + X2)x'>'^Y (6f) 



where 7 is the Euler gamma, and where we use v to 
mark contributions due to black-hole absorption effects; 
it is set to one to take those effects into account. 
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The TaylorT4 approximant is obtained by Taylor ex- 
panding the right hand side of Eq. (jSbp for dx/dt and 
truncating it at the appropriate reactive PN order. 
Therefore, to construct GW phase evolution in the Tay- 
lorT4 approximant that contains all the 3.5PN accurate 
non-spinning and 2.5PN accurate spin effects, the fol- 
lowing set of differential equations are numerically inte- 
grated: 



d(l){t) 



dt 

dx(t) dx{t) 



GM 



^3/2 



dt 



dx{t) I 
dt 



Ins 



dt 
16 c3 
5GM 
274229 



dxjt) 

NS H T — Is 



dt 



(7a) 
(7b) 



x^ i 1 



487 
168 



X + An X 



3/2 



254 



^5/2 



72576 21 
"178384023737 1712 



3353011200 
856 
105 



In (16 a;) 



105 
3310 
189 



1475 
192 

7/2 



(7c) 



dx{t) 
dt 


16 c3 r 47 
5GM \ 12 


Xl +X2 




+ 


I (Xi'+X2')4 


79 

-^XlX2 




+ 


1 . 

~32 + 


3x1=^ + 3 X2') ly 



8347 
2016 



(Xi +X2) 



.5/2 



(7d) 



The construction of the TaylorEt approximant requires 
PN accurate expressions for uji, in terms of £. The ra- 
diation reaction induced inspiral is incorporated by ex- 
pressing C in terms of £. The restricted PN waveform 
associated with the TaylorEt approximant reads 



h{t) cx £{t) cos2(/)(i) , 



(8) 



The PN accurate temporal evolutions for (j){t) and £{t), 
that includes all the 3.5PN accurate non-spinning and 
2.5PN accurate spin effects, are obtained by solving the 
following coupled differential equations. 



d(f>{t) 
dt 

^^b{t)\-NS 



UJbit) = Wh(t)|NS + Wfc(i)|s , 

12659 



^^3/2^ 37 
GM^ \ 32^ 



(^b{t)\s = 



GM 
205 
256 

,3 . r 7 



2048 



3016715 



__^3/2j i 

GM^ 1 4 



196608 
Xl 



(9b) 



X2 



(Xi +X2) 



^, 655 

£. 

^ 64 



t3/2 



Xl +X2 



^5/5 



(9c) 



dt 



dm ^ d^t) , 

dt 
16 

5GM 



NS 



dm, 



dt 



(9d) 



dm , 

dt 



Ins 



197 
336 



e + 47re 



3/2 



374615 2 299 
72576 ^ '^les'^^ 



5/2 



-I- 



3155 



dm^ 
dt 



1712 
T05" 
1712 
T05"^ 
16 

5GM 

7 f 2 
4 

~32 * 
9943 



In (4^^) 

e 



384 

4324127729 
82790400 
4155131 



96768 



7/2 



(9e) 



448 



(Xi 









Xl +X2 




111 


xl) + 


^X1X2 


+ X2 - 


^3x? + 3 


+ X2) 


^5/2 1^ 



t3/2 



(9f) 



where ^ = —2S/^<? and /i is the reduced mass. In this 
paper we keep d(j)/dt to its highest PN order and change 
PN orders of d^/dt to create various PN accurate Tay- 
lorEt approximants. The values of ^ corresponding to 
any initial and final GW frequencies can be numerically 
evaluated using the right-hand side of Eq. ([^ . This is 
possible due to the fact that for GWs from spinning com- 
pact binaries, having negligible orbital eccentricities, the 
frequency of the dominant harmonic is /gw = ^^6/2. 

We want to emphasize that the above mentioned Tay- 
lor approximants can only provide, from a strict PN point 
of view, fully 2.5PN accurate inspiral phase evolution for 
spinning black-hole binaries. This is because there are yet 
to be computed spin effects appearing at conservative and 
reactive 3PN and 3.5PN orders. Another point is that in 
these Taylor approximants, for the time being, we can 
only incorporate next-to-leading order spin-orbit effects 
and all other spin effects are mainly at the 'Newtonian or- 
der' (Newtonian order in the sense one refers to radiation 
reaction appearing for the first time at 2.5PN (absolute) 
order as Newtonian radiation reaction). Therefore, it is 
quite conceivable that smooth convergence to exact GW 
phase, reported in Ref. 48] for the TaylorEt approximant 
for non-spinning compact binaries, would not be present 
in its spinning counterpart. 



IV. NR-PN COMPARISON DURING INSPIRAL 

In this section we compare the phase and amplitude 
of the numerical waveforms with that predicted by post- 
Newtonian approximants. We have a lot of freedom in 
how to make such a comparison. Here we follow the pro- 
cedure we used in [l3| , which is also used in [3| : we line 
up the phase and frequency of the PN and NR waveforms 
when the gravitational- wave frequency AIlu is equal to 
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FIG. 2: An example of a phase comparison between NR and 
PN results. The case shown here is Si/M^ = 0.5, and the 
PN results were obtained with the TaylorTl approximant at 
2.5PN order. The accumulated phase error over the ten GW 
cycles before A4u = 0.1 is = —2.69 radians. 



0.1. This occurs several orbits before merger (the num- 
ber of orbits depends on the value of the black-hole spin, 
as suggested by Table at which time we expect post- 
Newtonian results to still be valid. We then compare the 
phase and amplitude at earlier times. When comparing 
results between different cases, we need to decide whether 
to compare over a period of time, a frequency range, or a 
number of GW cycles. The results do not change qualita- 
tively with different choices; we simply must make some 
choice; and the choice we make is to compare over a given 
number of GW cycles, which it seems to us will be most 
useful when later constructing hybrid waveforms. Our 
final procedure, then, is to compare the ten cycles before 
the GW frequency reaches Mw = 0.1. (This corresponds 
roughly to a black-hole orbital frequency of Mub — 0.05.) 



A. Phase 

As an example of how our comparison procedure 
works. Figure [2] shows the accumulated phase disagree- 
ment between the NR results for black holes with spins 
Si/Mf = 0.5, and the PN result using the Taylor Tl 
approximant at 2.5PN order. The phase disagreement 
is zero by construction when Mw = 0.1, and we choose 
i = at this point. As we progress backwards in time, 
the phase disagreement at first grows quadratically (since 
4'PN = 4'NR at t = 0, (A^) = 0, and so there cannot 
be any linear growth in A(/) at "early" times). A linear 
dependence soon develops, however, and the phase dis- 
agreement grows roughly linearly for the remainder of 
the comparison, back to t k, — 900M, when (f)NR has de- 
creased by 207r, or ten GW cycles. Note also that there 
are wiggles in the plot; these are due to the small eccen- 
tricity in the binary that was simulated numerically, and 
we expect that they introduce an overall uncertainty in 
the numerical phase of about 0.2 radians with respect to 
a binary with zero eccentricity (see Figure 15 in |13j). 

In the example shown in Figure [5] the accumulated 
phase error over the ten cycles before Mu = 0.1 is —2.69 



FIG. 3: The accumulated phase disagreement between NR 
and PN results over the ten cycles before Muj = 0.1. The five 
different spin values are Si/Mf = 0,0.25,0.5,0.75,0.85, and 
we compare with the three approximants TaylorTl, TaylorT4 
and TaylorEt. All approximants are calculated at 2.5PN or- 
der. 



radians. This is the most important piece of informa- 
tion we obtain from the plot, and rather than produce 
many very similar plots for different values of black- 
hole spin compared with different PN approximants at 
different PN orders, we summarize our results in Fig- 
ures [3] and m Here we see the accumulated phase agree- 
ment over ten cycles for five different values of spin, 
[Si/Mf = {0,0.25,0.5,0.75,0.85}), compared with three 
PN approximants, all calculated at 2.5PN and 3.5PN ac- 
curacy. The three approximants are the standard Taylor 
Tl approximant, the Taylor T4 approximant introduced 
in |14| . and the Et approximant introduced in [561 and 
compared in detail with the nonspinning case in [48| . The 
small eccentricity tends to increase (j^NRi so the points in 
Figures [3] and [4] are systematically too low by at most 0.2 
radians for the spinning cases, and we should also recall 
the numerical uncertainty in 4>mr of 0.25 radians. 

The main points we wish to emphasize in Figure [3] 
are that the accumulated phase disagreement between 
NR and 2.5PN TaylorTl results is roughly constant for 
all values of black-hole spin. This suggests that when 
producing hybrid waveforms for these cases, the same 
number of numerical cyles are needed for all spin values 
to produce hybrid waveforms of the same phase accuracy. 
On the other hand, the T4 and Et approximants have 
the advantage that the phase disagreement is less than 
it is for Tl, and decreases for higher spins, with the Et 
approximant performing best in the high-spin cases. 

Figure 2] shows a similar plot, but this time the ap- 
proximants are calculated up to 3.5PN accuracy. As 
pointed out in Section IIIIl the 3.5PN results are not 
3.5PN-accurate in all terms, and these results will change 
when all 3.5PN terms are known and included. 

We note once again that the phase disagreement be- 
tween the NR and TaylorTl phases is roughly constant 
for all values of spins. We also see that, as already seen 
in [l3 |. the TaylorT4 approximant agrees extremely well 
with the NR phase in the nonspinning case. However, 
this does not hold for the spinning cases, and for larger 
spins 3.5PN TaylorT4 performs worse than TaylorTl. 
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FIG. 4: The same comparison as in Figure |3] but this time 
the FN approximants are evaluated at 3.5PN order in those 
terms where this is possible (see text). 



The new TaylorEt approximant agrees extremely well in 
one case (when Si/Mf — 0.25), but gives a large dis- 
agreement for high spins. 

Our general conclusion from all of these comparisons 
is that, for the very small section of parameter space we 
have considered (equal-mass binaries with non-precessing 
spins in the "orbital hang-up" configuration), the phase 
agreement is roughly as good as in the nonspinning case; 
it certainly has not got dramatically worse, as one might 
have suspected. These numerical simulations can there- 
fore be used to produce hybrid waveforms of comparable 
phase accuracy to those that can be produced from our 
previous nonspinning data [1, [§, [l^ [sj] . 



B. Amplitude 

We now compare the amplitude of the {I — 2, to ~ 
2) mode of the NR with that calculated from the 
restricted PN approximation, which is the same as the 
amplitude from the quadrupole formula. The restricted 
PN amplitude of the GW strain h is given by Arcstrictod = 
(Mu;/2)^/'^/i?, where R is the distance from the source, 
and we have 



^ — ^rcstrictcd^ 



i<t>PN{t) 



(10) 



The phase (ppNit) depends on the PN approximant used 
in the previous section. However, if we differentiate h 
twice with respect to time to calculate and con- 

sider the amplitude as a function of GW frequency Mui, 
then we find that the choice of approximant has a neg- 
ligible effect on the resulting function. In other words, 
the GW strain amplitude Arostrictcd depends only on the 
GW frequency, and so does the restricted PN amplitude 
of r'I'4,22; to a good approximation. This point was also 



made in This fact allows us to make a comparison 



between the NR and restricted PN amplitudes that is in- 
dependent of the PN approximant that we used to calcu- 
late the phase — the results are independent of whether 
we use Taylor Tl, T4, or Et. 

In the nonspinning case we found that the NR and re- 
stricted PN amplitudes disagree by about (6 ± 2)% p^ . 



FIG. 5: The disagreement between restricted NR and PN 
amplitudes of r'I'4,22 as a function of GW frequency Mu). 
The lines below correspond, from bottom to top, to the cases 
Si/Mi = 0,0.25,0.5,0.75. (Every second line is dashed, to 
make them easier to distinguish.) In the nonspinning case the 
disagreement is roughly 6% over this frequency range. The 
disagreement increases as spin is added, and is about 11% for 
Si /Mi = 0.75. The large oscillations in this last case are due 
to the relatively high eccentricity of that system, e ~ 0.006. 



while [Tjl found 5% disagreement using numerical wave- 
forms with higher accuracy. Both sets of numerical re- 
sults found that the disagreement was roughly constant 
over the last ten cycles before Mu! = 0.1. 

The results for the spinning binaries we have consid- 
ered are shown in Figure [5l We see that for spinning 
binaries the amplitude disagreement increases. The plot 
shows results for spins Si/Mf = 0, 0.25, 0.5, 0.75. For the 
highest spin shown here, the amplitude disagreement be- 
tween NR and restricted PN is about 11%. In the case 
of Si/Alf — 0.85, fluctuations due to eccentricity make it 
difhcult to clearly measure the amplitude disagreement, 
and that case is not shown in the plot; however, we esti- 
mate the disagreement at around 12%, which is twice as 
high as in the nonspinning case. 

Higher-order PN amplitude corrections were able to 
improve the agreement in the nonspinning case to around 
2% using 2.5PN amplitude corrections [H, [l3|, and even 
better agreement was found when using 3PN amplitude 
corrections [l^ . It is clear from our current results that 
higher-oder PN amplitude corrections are even more cru- 
cial when spins are included. Higher-order corrections up 
to 2.5PN aheady exist [46, 54, 55], but 3PN amplitudes 
are not yet known for spinning binaries. These correc- 
tions may substantially improve the agreement with NR 
amplitudes, but we will postpone that study to future 
work. 



V. INSPIRAL AND MERGER 

We may ask how much the waveforms from spinning 
and nonspinning binaries differ. In data analysis appli- 
cations the overlap integrals between different waveforms 
determine how well a GW search will distinguish between 
them. If the phase of two waveforms can be lined up 
such that the waveforms are identical, then the overlap 
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between them is one. If the waveforms differ shghtly, the 
value of the overlap integral falls below one. In real-world 
data-analysis applications, the detector noise is included 
in the overlap integral, which leads to results that de- 
pend on the total mass of the binary. For the purposes 
of the purely qualitative illustration in this section, we 
do not use the detector noise, and we calculate the over- 
lap integrals in the time domain. Given some part of the 
waveform, we may calculate an overlap integral of the 
form, 

0{6t, S^) = -r^^ / ' hNs{t)h%{t + 6t, 54>)dt, (11) 
i\nsJ-^s Jti 

where "NS" indicates the waveform from a nonspinning 
binary, and "S" indicates one of the spinning-binary 
waveforms. The spinning waveform is time-shifted by 
some St, and phase-shifted by some Scj). The normaliza- 
tion factors are calculated by performing the integrals 

Nffs = r \hNs(t)\'dt, (12) 

Nl ^ r \hs{t)\''dt. (13) 
Jti 

For some choice of ti and t2, the time and phase shifts 
5t and Scj) can be optimized to find the largest possible 
overlap. 

Let us compare the nonspinning-binary waveform with 
that from a binary with Si/Mf — 0.75. We first choose 
a time interval that includes only inspiral, ti — 12QQM 
and t2 = 1700Af with respect to the nonspinning wave- 
form; this corresponds to roughly five GW cycles in the 
frequency range Mlo « 0.06 to Muj « 0.09. The maxi- 
mum overlap that can be achieved when comparing this 
small number of cycles with the S75 simulation is just 
under 0.99, meaning that the overlap is extremely good. 
This illustrates that during the inspiral, it is difficult to 
distinguish between the spinning and nonspinning wave- 
forms. Put another way, if one of these waveforms were 
detected, a large number of cycles would be required to 
estimate the black holes' spins. In practice, indeed many 
hundreds or thousands of cycles may be detected during 
the inspiral, and the accumulated phase difference be- 
tween spinning and nonspinning binaries may be more 
obvious. However, at present we must use PN approx- 
imants to model those hundreds or thousands of cycles, 
and it is not entirely clear how the variation in number 
of cycles when spin is added compares with the deviation 
of the PN approximants at different PN orders from the 
true solution. 

In addition, if a noticeable phase difference between 
spinning and nonspinning waveforms does accumulate 
over many hundreds of cycles, this might be partially 
compensated by rescaling the total mass of one of the 
binaries; a spinning binary may thus be detected as a 
nonspinning binary with the wrong mass. 

We now examine a time interval that includes the 
merger, ti = 1870M and t2 = 2020M. The maximum 



A 

/' \ 
/' ( 
/' 
1 

1 
















ji 
1 

11 
1 

— r 

— -S 


)12 

75 















1200 1300 1400 1500 1600 1700 
Time (M) 



0.4 
0.2 



-0.2 
-0.4 



\ 

\ 

\ 

I \ 
\ \ 


/' 1 

\ 1 


\ 

ll ' 
ll ' 
(1 ' 


ft 

1 h 


D12 

---- S75 


\ ' 1 
\ * / 
\ ^ /' 
\ \ 1 
\\ 1 

v// 




III 


ll 

Vi 


'i 1; 



1880 1900 1920 1940 1960 1980 2000 2020 
Time (M) 



FIG. 6: Spinning and nonspinning waveforms in the time in- 
tervals t = [1200, 1700] and t = [1870,2020], with the time 
and phase shifted so that they give the maximum value of 
the overlap integral (|ll|l . We can clearly see by eye that the 
overlap for the earlier time interval is much larger than that 
for the later interval. 



overlap that we can now achieve is less than 0.90: the 
waveforms differ far more in the merger phase. This is a 
far greater difference than we observe during the inspiral, 
and a far greater difference than we would expect to see 
even if many more inspiral cycles were available. 

Furthermore, if we were to adjust the total mass of one 
of the binaries, such that the overlap of the inspiral wave- 
forms was increased to almost one, this would not sig- 
nificantly improve the overlap of the merger waveforms. 
Therefore, detecting both the inspiral and merger may 
allow a significantly better estimate of the black holes' 
parameters. 

In Figure [6] we illustrate this point by showing the 
waveforms for the nonspinning and S'i/Mf = 0.75 simula- 
tions, over the two time intervals we just discussed, and 
lined up such that the overlap integral (|lip is a maxi- 
mum. It is immediately clear from the two plots that the 
overlap will be larger for the earlier time interval. 

We conclude from this simple analysis that, at least for 
binaries in an orbital hangup configuration, an accurate 
estimate of the black holes' spins will be far easier if the 
merger waveform is detected, rather than the inspiral. 
And, perhaps not surprisingly, the most accurate param- 
eter estimation would be possible if both the inspiral and 
merger were detected. A similar conclusion can be made 
when looking at higher modes, although all modes up to 
I = 4 have an amplitude no more than 10% that of the 
[l — 2,m = 2) mode in the cases we consider. 
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VI. DISCUSSION 

We have extended ourprevious comparison of NR and 
PN inspiral waveforms [13] to spinning binaries. In par- 
ticular, we considered equal-mass binaries whose black 
holes have equal spins oriented parallel to the orbital an- 
gular momentum. In these cases the spins do not precess 
and there is no net radiated linear momentum, but there 
is an "orbital hangup" effect that delays the merger as 
the spin is increased 15 1. 

We compared the accumulated phase difference be- 
tween NR and PN approximants over the ten cycles be- 
fore the GW frequency reaches Mcu = 0.1. We used the 
PN approximants TaylorTl, T4 and Et, at orders 2.5PN 
and 3.5PN. We conclude that, as in the nonspinning case, 
the growth of phase error is fairly small; for the number 
of cycles we considered, it is less than 2.5 radians for 
all 3.5PN approximants for spins up to Si/Mf — 0.85, 
and less than 1.5 radians for the TaylorTl approximant. 
The 3.5PN TaylorT4 approximant, which performed ex- 
tremely well in the equal-mass case [l3| is seen to perform 
much worse in the spinning cases. However, we stress 
that not all terms in the PN approximants are known 
to 3.5PN order. At 2.5PN order (the highest at which 
all terms are known), the Et approximant gives the best 
phase agreement. These results are summarized in Fig- 
ures [3] and [H 

We also compare the amplitude of the = 2,to = 2) 
mode of the NR r'i>4 with that calculated by the re- 
stricted PN approximation. In the equal-mass case the 
PN amplitude was larger than its NR counterpart by 



about 6%. In a binary whose black holes have spins 
Si/Mf = 0.85, the amplitude disagreement grows to 
around 12%. 

Finally, we gave an illustration of how the inspiral and 
merger waveforms differ between the spinning and non- 
spinning cases. Our simple analysis suggests that the 
spin is difficult to distinguish during inspiral, but has 
an extremely clear effect on the merger waveform. As 
such, we expect that estimating the spins of a binary's 
constituents should be much easier in cases where the 
merger is also detected. 

In the future we intend to extend this analysis to a 
larger sample of parameter space (in particular, more 
general spin configurations), and to explore the data- 
analysis implications of our results. 
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